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Abstract. In this paper, we examine the possibility of static, spherically symmetric 
. gravitational geons on a 3 dimensional brane embedded in a 4+1 dimensional space- 

time. We choose a specific gu for the brane-world space-time metric. We then calculate 
g rr analytically in the weak field limit and numerically for stronger fields. We show 
that the induced field equations on the brane do admit gravitational geon solutions. 
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1. Introduction 

> 

in 

\q \ General Relativity in 3+1 dimensions does not admit static, nonsingular, asymptotically 
flat and topolotically trivial vacuum space-time solutions [IHH]. In this article, we 
investigate whether such configurations of the gravitation field, or gravitational geons 
[IHH], are possible on the brane in the Randall and Sundrum model [7J[8] in which the 
three spatial dimensions that we experience are a 3 dimensional brane embedded in a 
4+1 dimensional space-time. The effective 3+1 dimensional Einstein field equations on 
the brane for the Randall and Sundrum model were derived by Shiromizu, Maeda, and 
H ' Sasaki M. 

Gravity on the brane is prevented from leaking into the extra dimension at low 
energies by a negative bulk cosmological constant [10J. This negative cosmological 
constant can be offset by a positive brane tension such that there is a zero cosmological 
constant on the brane. In the case of a gravitational geon solution in which no matter 
is present, there would nevertheless be a generally non-zero effective energy-momentum 
tensor induced on the brane corresponding to 

( %, = ^C^VaV'KK (i) 

where ^S^ v is a 4-dimensional projection of the 5-dimensional Weyl tensor ^C^ pa , n a 
is normal to the brane and = g pv — r]^r] u is the induced metric on the brane [9|[TT]. 

In their investigation of black holes on the brane, Dadhich, Maartens, Papadopoulos 
and Rezania showed that the Reissner-Nordstrom metric is an exact solution to the 
effective 3+1 dimensional Einstein field equations on the brane and can be interpreted 
as a black hole without electric charge but with a tidal charge that arises from the 
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solution to the Einstein field equations in 4+1 dimensions [12J. They further showed 
that, in the induced metric for their exact black hole solution, the free gravitational 
field in the bulk contributes a negative effective energy density on the brane and tends 
to strengthen the gravitational field. 

In this paper, we investigate whether the free gravitational field in the bulk modifies 
the gravitational field on the brane so that gravitational geons are possible on the brane. 
Our approach is to look at a specific static, spherically symmetric solution of the form 

ds 2 = -B(r)dt 2 + A(r)dr 2 + r 2 (d6 2 + sin 2 Odcf 2 ) . (2) 

For this metric, the 4 dimensional Ricci scalar on the brane is given by [11 J 



B" B' (A' B'\ 2 (B' A'\ 2 2 



AB 2AB \ A B Ar \ B A Ar 2 r 2 



where the prime indicates differentiation with respect to r. 

Our approach is to choose B{r) and solve R = in <^ for A{r). Substituting a 
solution to R = into the Einstein field equations for 3+1 dimensional general relativity 
gives a trace free energy-momentum tensor 

that satisfies 

V%, = . (5) 

For a vacuum solution on the brane, the projection of the 5-dimensional Weyl tensor 
([U) satisfies 



= , (6) 

allowing us to make the correspondence SnT^ = —S^ and treat R = along with 
equations (j4]) and ([5]) as a closed system of equations on the brane. The effect is that 
(from [T2]) 

"a stationary general relativity solution with trace-free energy-momentum 
tensor gives rise to a vacuum brane-world solution in 5-dimensional gravity. " 

Examples of other papers that use this correspondence to investigate vacuum solutions 
on the brane with R = (or R = 4A if one assumes a non-zero cosmological constant 
on the brane) include [!Tl[T4Tfl8] . 

In this paper we will show that gravitational geons are possible on the brane. 

We are not interested in space-time solutions with singularities so we examined the 
Kretschmann scalar for divergent behaviour in order to identify candidate space-times. 
The Kretschmann scalar is given by 

R„v P aR^ pa = + 8Kl + 8K 2 + AK 2 A (7) 
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where 

I (B» (B'f A'B'\ 
Ki -a{2B~^~AAb) ' (8) 

k ^4b^ (9) 

K 3 = ^~ r and (10) 
From ( flOl) and ( fTTj) we must have limA(r) = 1 + 0(r n ) with n > 2 to avoid 

7 >0 

5'(r) 

divergence of the Kretschmann scalar. Similarly, from (T9l) it follows that lim — — — = 

r-»o B(r) 

0(r m ) with m > 1. 

In this paper, we take 

B( r ) = 1 - Q 2mr2 „ , (12) 

for which B'(0) = 0, -B(O) = 1, both B(r) and B'{r) are continuous for all r > and 
B(r) is Schwarzschild in the limit as r — > oo. We also consider general functions B(r) 
in the weak field limit. 

The solutions that we investigate have A(r) and B(r) greater than zero and finite 
for all r > 0. This ensures that the spacetime will be nonsingular and without horizons. 

2. Outline of Our Approach 

For static solutions, we have R = which can be written as 

_ {2BB"r 2 - [B'fr 2 + ABB'r + AB 2 )A - AB 2 A 2 
~ Br(B'r + AB) ' ^ ' 

We did not find a general solution to f TTS]) for our choice of B{r). In the absence of 

a general solution, our approach is as follows 

• Solve f fT3|) in the weak field limit. 

• Use directional field plots (generated using Maple) and iterative numerical analysis 
to discover the behaviour of A(r). 

• Find solutions near potential singular points (i.e. at r = and B'r + 4B = 0) in 
ffT3"|) using series approximation. Specifically, we write (fT3~j) in the form 

A\r) = f(r)A(r) + g(r)A 2 (r) . (14) 

Letting F'(r) = f(r) then ( Tl4"|) has the solution [H] 

A ^ = e-nr) J g(r) e nr)dr ' (15) 

• Examine the behaviour of the Kretschmann scalar R^ pu R^ up(T at potential singular 
points. 
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3. The Weak Field Approximation 

In the weak field, taking B(r) = 1 + b(r) and A(r) = 1 + a(r) with \b(r)\ << 1 and 
\a(r)\ « 1, ffTBT) can be written as 

, b"r 2 + 2b'r-2a 
a fa . (16) 

This has the solution (where c is a constant of integration) 

a = T + ; . (17) 

Noting that c = to ensure that A is finite at r = 0, any choice of B(r) for which b'r is 
small for all r would give suitable geon solutions in the weak field. 
For our choice of B(r) this becomes 

mr 5 — 4m 2 l 2 r 2 

a= (r 3 + 2m/ 2 ) 2 ' (18) 

For r << (ml 2 ) 3 we have B fa 1 + (jj and A ~ 1 — (y) . A(r) therefore satisfies 
the requirements that A(0) = 1 and limA(r) = 1 + 0(r n ) with n > 2. It is worth 

r— >0 

noting that for large r, A differs from Schwarzschild with A fa 1 + — . Furthermore, the 
conditions for the weak field, 



\a[r)\ 



mr 



5 



4m 2 l 2 r 2 



(r 3 + 2m/ 2 ) 



2 



« 1 (20) 



are both met for j << 1 and for y « 1. The solutions (see figH]) correspond to 
gravitational geons. 



4. General Behaviour 

We did not find a general solution to ( fl3l) for strong gravitational fields. We examined 
the behaviour of A(r) using a fourth order Runge-Kutta numerical iteration. Using this 
method, we were unable to begin at exactly A(0) = 1. However, A fa 1 — (j\ for 
small r. Thus, if we choose an iteration length Ar so that = 10~ 5 then A(Ar) fa 1 
to within 10~ 10 (independent of m). Also, having made the observation from direction 
field plots that the solution is insensitive to small changes in the initial value of A, we 
chose Ar appropriately and used A(Ar) = 1 as the starting point of the iteration. 

We are interested in solutions for which both A(r) and B(r) are greater than zero 
and finite for all r > 0. The zeros of B can be determined by noting that B(0) = 1, 
lim B(r) = 1 and that B' = at r = (Ami 2 ) 5 where B has its minimum value of 

1 _ |(4f )§. Thus 
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no zeros 


if 


m 
I 


< 


V27 
4 




one zero 


if 


in 
I 




4 


(21) 


two zeros 


if 


m 
I 


> 


V27 
4 ' 





We proceeded to examine the numerical iteration of A(r). We found that it works 
well provided we choose m and I for which B'r + 45 ^ for all r > but that it is 
unable to navigate instabilities in the solution otherwise. The zeros of B'r + 45 are 
found by observing that B'r + 4B\ r=0 = 4, lim B'r + AB = 4 and that J^(5'r + 45) = 

at r = (2ml 2 ) 3 where B'r + 45 has its minimum value of 4 — (^p 1 )^. Thus 



no zeros if j < p 
5V + 45 has { one zero if f = p (22) 
two zeros if f > p . 

Noting that there will be choices of m and Z for which there will be zeros of B'r + 45 
but no zeros of B, we investigated the behaviour of A(r) around the zeros of B'r + 45. 
In general, from (1T4T) . we have 

, 255"r 2 - (5') 2 r 2 + 455'r + 45 2 , , s 

/(r) = 5r(5V + 45) ^ (23) 

-45 , , 

r[B'r + 45) 

Taking x = r — r where B'r + 45| r=ro = we get series approximations for ( 1231) and 

/fx) = h «2 + a^x + 0(x 2 ) and (25) 

x 

0(ar) = ^ + f3 2 + fax + 0(x 2 ) . (26) 
x 



A i x ) = -m /3 ia2 - feQl , m — , o , i — ( A " const - of integration). (27) 



In this case, for ot\ £ {0, —1, —2, ...}, solutions to (JT5]) are of the form 

-1 

ai ai(oi+l) 1 I I V ' II 

For ai G {0, — 1, —2, ...} we get a different solution for each choice of a\. For example, 
choosing a.\ — — 1 gives 

A fx) = — - 77: r : 1 1 —, — 1 (A - const, of integration). (28) 

w -f3 l + (f3 2 + f3 l a 2 )xln\x\+(3 1 a 2 x + \x + 0(x 2 ln\x\) K ' K ' 

In each case, the dominant term in the denominator as x — > results in the same 
outcome as (127]) that 
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— - if ai < or A = 
\imA(x) = l &\ (29) 

x_> ° ( if ai > and A ^ . 

We calculated values for a± and fi\ for each zero of B'r + AB given a variety of 
values of m and / (see figs. [2] and [3]). Having verified our results for a range of choices 
0.02 < I < 750, we found that at each zero, the values of a\ and (3i depend on — but 
not on m and I individually. Knowing that zeros of A(r) correspond to singularities in 
the space-time, we sought to determine where a% = in terms of y. 

With x = r — To, we have B'r + AB x[B"r + 5B') \ r=ro in the vicinity of each zero 
of B'r + AB so 



a i 



2B "r-^ + AB' + ^ 

B r 



B"r + 5B' 
2B"r + IB' 



=''o 



B"r + 5B' 



For «i = 0, substituting in ( fl2l) . we get the following quadratic in r 3 



r D + Uml 2 r 3 - 24m z l 



,2,4 







(30) 



(31) 



with one positive real root corresponding to r ~ (1.544m/ 2 )s. Substituting this into 
£?V + AB = yields j ~ 1.202. Based on this (and again looking at figs. [2] and [3]), we 
know that for solutions with y > p, there will be two zeros of B'r + AB, that the first 
will have ct\ > for y > 1.202 and the second will always have a\ > 0. 

We know that A' — >■ ±oo corresponds to singular behaviour in the space-time. 
Differentiating (1271) gives 



A' (a;) 



filQ2-fe«l 



«i A|x| 



+ 



Pi 
ai 



/3ia 2 -feQ] 
Qi(ai + 1) 



X\x\ 



+ 



(32) 



lim A'(x) 



/3i 2 («i + 1) 
— s^n(aiAx)oo 




if ai < — 1 or A = 

if —1 < Qfi < 1 and A ^ 
if «i > 1 and A 7^ . 



(33) 



We know from ( 1291) that we need not consider «i > (unless A = 0), but the 
behaviour of A(r) does depend on whether ot\ is greater or less than -1 at the zeros 

—1 in terms of and 



B'r + AB. We therefore sought to determine where ct\ 
numerically found this corresponds to y ~ 1.191. 

Now we found that at each zero, the values of ct\ and (3i depend on y but not 
on m and / individually. That is not true of the values of a 2 and (3 2 . However, we did 
determine (again verifying our results for a range of choices 0.02 < / < 750) that sgn(a 2 ), 
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sgnlfa) and the ratio ^ at each zero depends on y but not on m and Z individually. 

We can therefore say that the sign of ^1^+1) = — /3^(«i+i) — depends on y but 
not on m and I individually. 

Having determined how the choice of y affects the behaviour of A and A' at each 
zero of B'r + AB, we looked at each distinct range of y values considering both A = 
and A 7^ solutions. For each case, we sought to determine if our solution corresponds 
to gravitational geons. 

4-1- y < p - Solutions correspond to gravitational geons 

For these values of y, B'r + AB has no zeros. Consequently, all solutions with A(0) = 1 
correspond to gravitational geons. Choosing I = 1 we confirmed the behavior of A using 
iterative numerical analysis (see figJl]). Also evident is the progression from the weak 
field solution for y << 1 (figHJ) to the piecewise constructed solution for p < — < 1.191 
(figJHJ) which we discuss in the following section. 

4.2. p < y < 1.191 with A = at the second zero of B'r + AB 

At the first zero of B'r + AB, A = > and A' = ai{ t^~fu 2) > for values of f 
in this range. There is therefore no singular behaviour and no horizon indicated at the 
corresponding spherical hypersurface. 

At the second zero of B'r + AB, A = indicating singular behaviour in the space- 
time except when A = 0. In that case, A = > and A' = ^g^yf^ > at the 
second zero. Consequently, no singular behaviour or horizon is indicated. 

Whereas our analysis of the Kretshmann scalar shows the possibility of solutions 
corresponding to gravitational geons with A = at the second zero, instability at 
each zero prevents us from confirming any particular solution using iterative methods 

— Oil 

(see fiej6l). However, we have been able to confirm that lim A = and lim A' = 

x w x->0 Pi x^O 

1 1 — — ^ 2 around each zero for each section constructed numerically (see figs. [7J 
Pi{atx + 1) 

mini). 

We piecewise constructed A around each zero of B'r + AB using the numerically 
constructed sections (see figJSj . The resulting solution is sewn together at the spherical 
hypersurfaces corresponding to each zero. 

The sewing together of two manifolds may induce a surface energy-momentum 
tensor [HIED] 

S tlv = ^{[K llv ]-[K]h llv ) (34) 

on the surface where they are joined. Here K^ u is the extrinsic curvature of the surface, 
K = M , h^ v is the induced metric on the surface and [K^] denotes the jump in K^ v 
across the surface. 
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Now, K^y depends on the metric and its first derivatives [2T]. Given that 
A, B and B' are all continuous in our piecewise constructed solution, there will be a 
non- vanishing only if A' is discontinuous across either hypersuface corresponding to 
a zero of B'r + AB. Our analysis shows that A' is continuous implying that S^ u = 0. 

It should be noted that there could be contributions to S^ v from the Weyl term 
([T]). If that were the case, S^ v = in our piecewise constructed solution would imply 
that there is a localized gravitational source (such as matter) at the hypersurface and 
the stresses contributed by each cancel one another exactly. Since the space-time off the 
brane is not known, it is not possible to check to see if such contributions exist. 

If there are no contributions to S^ v from the Weyl tensor, the piecewise constructed 
solutions correspond to gravitational geons. 

4-3. Choices of y for which no gravitational geon solutions exist 

Our analysis shows that the Kretshmann scalar diverges at the spherical hypersurfaces 
corresponding to the zeros of B'r + AB for values of y > 1.202. We have therefore 
determined that solutions that correspond to gravitational geons do not exist for these 
values of y. 

In the range 1.191 < y < 1.202, divergence of the Kretshmann scalar rules out 
the possibility of gravitational geon solutions, except in the special case that A = at 
both zeros of B'r + AB. However, direction field plots indicate that these solutions are 
inconsistent with our boundary conditions. 

As discussed in the previous section, gravitational geon solutions are again ruled 
out by divergence of the Kretshmann scalar in the range |f < y < 1.191, except in the 
special case that A = at the second zero of B'r + AB. 

There are therefore no solutions that correspond to gravitational geons possible for 
y > §f with the exception of A = at the second zero of B'r + AB for p < y < 1.191. 

5. Conclusion 

In this paper, we investigated a particular static, spherically symmetric 3+1 dimensional 
space-time for behaviour consistent with that of gravitational geons on a brane 
embedded in a 4+1 dimensional space-time. In the absence of a general solution, we 
investigated the behaviour of the space-time using a variety of methods and were able to 
ascertain that, for a particular set of parameter choices, gravitational geons are possible. 

The significance of the result presented in this paper is that, in general, for brane- 
world space-times gravitational geons are possible. We investigated solutions in the 
weak field limit and showed that gravitational geons will exist as long as B'r is small for 
all r > 0. We also investigated the specific function B{r) = 1 — ^^mi 2 anc ^ snowe d that 
gravitational geon solutions exist for all parameter choices such that y < §f • These 
solutions include both weak field and strong field geons. 

If our universe corresponds to the Randall- Sundrum model, then gravitational 
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geons could represent a form of dark matter (this possibility has also been discussed 
by Sones [6] for quantum geons with a Klein-Gordon field). 
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Figure 1. B(r) and A(r) plotted against radius r for I = 1, m = 0.01. A 4th 
order numerical iteration of (|T3j) is labelled A(r) and the analytical weak field result 
is labelled weakA(r). These two approaches give almost identical results for y << 1. 
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m/l 



Figure 2. ai and (3i plotted against y for the first zero of B'r + AB. 
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Figure 3. a\ and (3i plotted against y for the second zero of B'r + 4B. 
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Figure 4. Function A(r) plotted against radius r. Here we choose I = 1 and chart 
a variety of values for m in the specified region. A 4th order Runge-Kutta numerical 
iteration is used to generate each result. 




Figure 5. The function A(r) plotted against radius r for I = 1, m = 1.189. A 4th 
order Runge-Kutta numerical iteration was unable to navigate the zeros of B'r + 4_B 
(at r — 1.281, 1.392), so the result is constructed piecewise between these points. 
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Figure 6. This direction field plot of A(r) against radius r shows that A(r) is 
numerically unstable at the zeros of B'r + AB. Here we choose I = 1, m = 1.189 
such that the zeros are found at r = 1.281, 1.392. 
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Figure 7. Forward iteration of A(r) plotted from r = to the first zero of B'r + 4£> 
for I = 1, m = 1.189. A — > ^-p- w 1.38 as the solution approaches the first zero. 




Figure 8. Reverse iteration of A(r) plotted from the second zero of B'r + AB to the 
first for / = 1, m = 1.189. 
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Figure 9. Forward iteration of A(r) plotted from the second zero of B'r + 4i? for 
l = l,m= 1.189. 



